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p ^. Abstract 

I A closed formula is obtained for the infrared singularities of dimensionally regularized, 

massless gauge-theory scattering amplitudes with an arbitrary number of legs and loops. 
, It follows from an all-order conjecture for the anomalous-dimension matrix of n-jet oper- 

\ ators in soft-collinear effective theory. We show that the form of this anomalous dimen- 

■ sion is severely constrained by soft-collinear factorization, non-abelian exponentiation, 

,—1 . and the behavior of amplitudes in collinear limits. Using a diagrammatic analysis, we 

'""J I demonstrate that these constraints imply that to three-loop order the anomalous dimen- 

' sion involves only two-parton correlations, with the possible exception of a single color 

. structure multiplying a function of conformal cross ratios depending on the momenta of 

I four external partons, which would have to vanish in all two-particle collinear limits. We 

' suggest that such a function does not appear at three-loop order, and that the same is 

true in higher orders. Our formula predicts Casimir scaling of the cusp anomalous di- 
mension to all orders in perturbation theory, and we explicitly check that the constraints 
I exclude the appearance of higher Casimir invariants at four loops. Using known results 

for the quark and gluon form factors, we derive the three- loop coefficients of the l/e"' 
pole terms (with n = 1, . . . , 6) for an arbitrary n-parton scattering amplitude in massless 
QCD. This generalizes Catani's two-loop formula proposed in 1998. 
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1 Introduction 



The origin and structure of ultraviolet (UV) divergences in quantum field theories is well 
understood. They can be absorbed into a renormalization of the parameters of the theory. 
The fact that physical results must be independent of the UV regulator introduced in in- 
termediate steps of a calculation gives rise to powerful constraints, which are summarized 
by the renormalization-group (RG) equations of the theory. Perturbative results for on-shell 
scattering amplitudes in theories with massless fields also contain infrared (IR) singularities, 
which originate from loop-momentum configurations where particle momenta become soft or 
coUinear. These singularities cancel in physical observables, which also include real radiation 
and are insensitive to soft and coUinear emissions [1, 2]. 

In a recent letter [3], we have shown that the IR singularities of on-shell, n-particle scat- 
tering amplitudes in massless QCD are in one-to-one correspondence with UV divergences of 
operators defined in soft-coUinear effective theory (SCET) [4-7]. This implies that they can 
be analyzed using standard methods of operator renormalization. In particular, the IR diver- 
gences of n-point scattering amplitudes can be absorbed into a multiplicative renormalization 
factor Z and are fully determined by an anomalous dimension F. Both the Z'-factor and the 
anomalous dimension are matrices in color space, i.e., they mix amplitudes with the same 
particle content but different color structure. In dimensional rcgularization, IR singularities 
manifest themselves as poles in e = (2 — rf/2). The renormalization conditions imply that the 
higher 1/e" poles are determined by the single-pole coefficient, from which one obtains the 
anomalous dimension. In contrast to traditional effective theories, SCET operators are nonlo- 
cal along the direction of large momentum flow. This induces a dependence of the anomalous 
dimension V on the momentum transfers and translates into the presence of Sudakov double 
logarithms in high-energy scattering processes. 

The structure of the effective theory imposes non-trivial constraints on the form of the 
anomalous dimension F. SCET contains a separate set of coUinear fields for each direction 
of large energy fiow, and a single set of soft fields mediating interactions among the different 
coUinear sectors. Up to contributions suppressed by powers of the large momentum transfers, 
the different coUinear fields interact only through soft gluon exchange. The matrix elements 
of n-parton operators factor into jet functions associated with each set of collincar fields and a 
soft function which arises from interactions among the different jets. Since different coUinear 
fields do not interact, the jet functions are color diagonal. Non-trivial momentum dependence 
and color correlations of the anomalous dimension arise from soft interactions. The fact that 
the soft and jet functions must conspire to produce an anomalous dimension depending only 
on hard momentum transfers leads to strong constraints on the form of the allowed terms. 

Due to the eikonal form of soft interactions, they can be represented as Wilson lines along 
the directions of the external colored particles. This leads to further constraints implied by 
the non-abclian exponentiation theorem [8, 9]. In QED, eikonal identities imply that higher- 
order soft radiation is obtained by exponentiating the leading-order contribution [10]. In 
QCD this is no longer true, but only a particular set of color structures can contribiitc to the 
soft anomalous dimension. We argue that their coefficients are related to each other by the 
structure of eikonal interactions. Additional constraints are obtained by considering limits in 
which some of the partons involved in a scattering process become coUinear. In the simplest 
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case where two external momenta become coUinear, an n-parton amplitude reduces to a sum 
of {n — l)-parton amplitudes multiplied by splitting amplitudes. These splitting amplitudes 
and their anomalous dimensions can only depend on the momenta and quantum numbers of 
the particles involved in the splitting process. This leads to a constraint on the difference of 
the anomalous dimensions of the amplitudes before and after the splitting. 

In this paper we derive the most general form of the anomalous dimension T compatible 
with the above constraints. We argue that to all orders in perturbation theory the anomalous 
dimension has the same structure as at one-loop order, featuring only two-particle correlations 
in colors and momenta. Our result is semi-classical in that it describes color-dipole interactions, 
which probe the momenta and color charges of pairs of external particles. This is reminiscent 
of Low's theorem, which relates the cross section for soft photon emission to the classical 
cross section times a factor depending on the electric charges and momenta of the external 
particles [11]. In our case, the quantum structure of the underlying field theory shows itself via 
the coefficient of the color-dipole operator, which we relate to the universal cusp anomalous 
dimension of light-like Wilson loops, and via anomalous dimensions for quark and gluon fields. 
In practice, these anomalous dimensions can be determined at three-loop order using existing 
results for the quark and gluon form factors in QCD. However, our analysis is completely 
general and extends to arbitrary gauge theories with massless fields. 

Starting at three-loop order, the analysis of the structure of the anomalous dimension 
becomes non-trivial. Based on a general study of the constraints following from soft-coUinear 
factorization and the non-abelian exponentiation theorem, we show that in the most general 
case two additional color structures beyond those predicted by our simple formula could appear 
at three-loop order. One of them describes interactions among four different partons. We also 
show that at four-loop order a new structure could arise, which would violate Casimir scaling of 
the cusp anomalous dimension. However, imposing the correct behavior of n-parton scattering 
amplitudes in the limit where two partons become coUinear eliminates all of these additional 
structures. We also present a second argument for the absence of the additional terms, based 
on the simple form of color-symmetrized soft-gluon interactions, which suggests that these 
new structures vanish due to color conservation. In this way our formula is established at 
three-loop order. We note, however, that our arguments could be circumvented by functions 
of conformal cross ratios of four parton momenta that vanish whenever two partons become 
coUinear. At three-loop order a single color structure involving such a function is allowed by 
non-abelian exponentiation. For our conjecture to be valid, this function must vanish. The 
observation of Casimir scaling of the four-loop cusp anomalous dimension is unaffected by this 
caveat. 

The arguments presented in this work suggest that our conjecture may hold to all orders 
in perturbation theory. This would imply a new set of exact relations among amplitudes in 
perturbativc quantum field theory and will hopefully shed new light on their deeper structure. 
The knowledge of the structure of IR singularities, and of the associated anomalous- dimension 
matrix, is also of significant practical interest for collider physics. It is a necessary ingredient 
for the resummation of large logarithms in jet-production processes. Solving RG equations in 
SCET resums perturbativc logarithms of ratios of the invariant masses of the produced jets to 
the large momentum transfers involved in their production. As an application of our results, 
we predict the complete structure of IR singularities for n-parton scattering amplitudes in 
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massless QCD at the three-loop order. In recent work, the four-gluon amphtude in J\f — A 
super- Yang-Mills theory (SYM) was studied at three-loop order and expressed in terms of a 
small set of master integrals [12]. Once these integrals have been evaluated analytically or 
in numerical form, this calculation will provide a stringent test of our predictions. Another 
interesting implication of our analysis is the prediction that the cusp anomalous dimensions 
of quarks and gluons should obey Casimir scaling to all orders in perturbation theory, i.e., 
they should be equal to the quadratic Casimir operator Cr in the fundamental or adjoint 
representation times a universal coefficient. While Casimir scaling has been shown to hold up 
to three-loop order by explicit calculation [13], this prediction is highly non-trivial in view of 
the expectation that this scaling should no longer hold non-perturbatively, at least not for the 
finite parts of Wilson-loop expectation values. Already a long time ago, Prenkel and Taylor 
argued that Casimir scaling would be inconsistent with expectations about the area law for 
matrix elements of Wilson loops giving rise to confinement [9] . More recently, investigations 
of high-spin operators in string theory using the AdS/CFT correspondence [14-16] have found 
a strong-coupling behavior that is inconsistent with Casimir scaling [17-19]. 

We begin in Section 2 by recaUing the connection between on-shell scattering amphtudes 
and Wilson coefficient functions in SCET [3]. We then solve the RG equation for the renormal- 
ization factor Z in terms of an integral over the anomalous- dimension matrix of n-jet SCET 
operators. We derive the three-loop expression for Z in terms of known anomalous- dimension 
coefficients and show that our approach reproduces Catani's result for the divergences of 
two-loop amplitudes [20]. In Section 3 we discuss the structure of SCET for n-jet processes, 
explicitly construct the necessary operators and show that their matrix elements factor into jet 
and soft functions. We stress the importance of soft operators built out of n light-like Wilson 
lines, which describe the color and momentum correlations in the anomalous- dimension ma- 
trix. A detailed discussion of the arguments supporting our conjecture for the structure of the 
anomalous-dimension matrix is presented in Section 4. We recall important facts about the 
renormalization of Wilson loops with cusps and cross points, the non-abelian exponentiation 
theorem, and soft-collinear factorization in SCET. In Section 5 we study the implications of 
the known behavior of scattering amplitudes in the limit where two of the external partons 
become coUinear. We find that the n-jet anomalous-dimension matrix can be decomposed in 
this limit into the sum of an {n — l)-jet anomalous dimension and the anomalous dimension 
of the splitting amplitudes, whose all-order form we derive. The nature of this decomposition 
imposes another strong constraint on the momentum and color structures that can appear in 
the anomalous dimension. In Section 6 we perform an explicit analysis of the structure of the 
soft anomalous-dimension matrix to three-loop order. We first list all structures allowed by 
non-abelian exponentiation and then impose the constraints from soft-collinear factorization 
and from two-particle coUinear limits. These constraints eliminate all additional terms with 
the exception of a single color structure, which is of subleading order in the ^ oo limit 
and compatible with the constraints if it is multiplied by a function of conformal ratios which 
vanishes in all coUinear hmits. We also show that the constraints enforce Casimir scaling 
of the cusp anomalous dimension to four loops. Our conclusions are presented in Section 7. 
Perturbative results required to cvahiatc our formulae at three-loop order and some comments 
on the behavior of our color structures in the large- A^c limit are compiled in two appendices. 
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2 IR factorization and RG invar iance 



The key observation of our letter [3] was that the IR singularities of on-shell amplitudes in 
massless QCD are in one-to-one correspondence to the UV poles of operator matrix elements 
in SCET. These poles can therefore be subtracted by means of a multiplicative renormahza- 
tion factor which is a matrix in color space. Specifically, we have shown that the finite 
remainders of the scattering amplitudes can be obtained from the IR divergent, dimensionally 
regularized amplitudes via the relation 

|A^n(M,/^)) = lini Z-^(e,M,//) |Mn(e,M)) . (1) 

Here {p\ = {pi, . . . ,j9„} represents the set of the momentum vectors of the n partons, and 
/I denotes the factorization scale. The quantity |A^„(e, {p})) on the right-hand side is a 
UV-renormalized, on-shell n-parton scattering amplitude with IR singularities regularized in 
d — A — 2e dimensions. After coupling constant renormalization, these amplitudes are UV 
finite. Apart from trivial spinor factors and polarization vectors for the external particles, the 
minimally subtracted scattering amplitudes \-M.n{{p}, f^)) on the left-hand side of (1) coincide 
with Wilson coefficients of n-jet operators in SCET [3], to be defined later: 

|A4„({p}, //)) = X [on-shell spinors and polarization vectors] . (2) 

We postpone a more detailed discussion of the effective theory to Section 3 and proceed to 
study the implications of this observation. 

To analyze the general case of an arbitrary n-parton amplitude, it is convenient to use the 
color-space formalism of [21, 22], in which amplitudes are treated as n-dimensional vectors 
in color space. Tj is the color generator associated with the i-th parton in the scattering 
amplitude, which acts as an SU (iVc) matrix on the color indices of that parton. Specifically, 
one assigns {T°-)ap = f^p for a final-state quark or initial-state anti-quark, {T^"')af3 = —tpa 
a final-state anti-quark or initial-state quark, and {T^^)bc — —if"'^'^ for a gluon. We also use 
the notation Ti ■ Tj = T"- T"- summed over a. Generators associated with different particles 
trivially commute, Tj ■ Tj = Tj ■ Tj for i ^ j, while T^ = Ci is given in terms of the quadratic 
Casimir operator of the corresponding color representation, i.e., Cg = Cq = Cp for quarks or 
anti-quarks and Cg — Ca for gluons. Because they conserve color, the scattering amplitudes 
fulfill the relation 

~ Il°|X„(e,{p})> = 0. (3) 
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It follows from (1) that the minimally subtracted scattering amplitudes satisfy the RG 
equation 

-f- \Mn{{p}, /.)) = r({p}, /.) \Mn{{p}, 1^)) , (4) 

flln/x - - - 

where the anomalous dimension is related to the Z-factor by 

r(M, lA = -z-\e, M, ii) ^ z(6, M, //) . (5) 
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The formal solution to this equation can be written in the form 



Z{e,{p},^J) = Pexp 



(6) 



where the path-ordering symbol P means that matrices are ordered from left to right according 
to decreasing values of //'. The upper integration value follows from asymptotic freedom and 
the fact that Z = 1 + 

In Section 4, we will discuss theoretical arguments supporting an all-order conjecture for 
the anomalous-dimension matrix presented in [3], which states that it has the simple form 



where Sij = 2aijPi ■ pj + iO, and the sign factor aij = +1 if the momenta pi and pj are both 
incoming or outgoing, and aij = — 1 otherwise. Here and below the sums run over the n 
external partons. The notation (ii, ik) refers to unordered tuples of distinct parton indices. 
Our result features only pairwise correlations among the color charges and momenta of different 
partons. These are the familiar color-dipole correlations arising already at one-loop order from 
a single soft ghion exchange. The fact that higher-order quantum effects do not induce more 
complicated structures and multi-particle correlations indicates a semi-classical origin of IR 
singularities. Besides wave-function-renormalization-type subtractions accomplished by the 
single-particle terms 7*, the only quantum aspect appearing in (7) is a universal anomalous- 
dimension function 7cusp related to the cusp anomalous dimension of Wilson loops with light- 
like segments [23-25] . The three anomalous- dimension functions entering our result are defined 
by relation (7). They can be extracted from the known IR divergences of the on-shell quark 
and gluon form factors, which have been calculated to three-loop order [26-28]. The explicit 
three-loop expressions are given in Appendix A. 

Concerning the form of (7), we note that a conjecture that an analogous expression for 
the soft anomalous-dimension matrix (sec Section 4.4 below) might hold to all orders was 
mentioned in passing in the introduction of [12], without presenting any supporting arguments. 
In a very recent paper, Gardi and Magnea have analyzed the soft anomalous- dimension matrix 
in more detail and found that (7) is the simplest solution to a set of constraints they have 
derived [29]. However, they concluded that the most general solution could be considerably 
more complicated. Indeed, we emphasize that as a consequence of our result some amazing 
cancellations must occur in multi-loop calculations of scattering amplitudes. At L-loop order 
Feynman diagrams can involve up to 2L parton legs, while the most non-trivial graphs without 
sub divergences can still connect (L -|- 1) partons. We predict that these complicated diagrams 
can be decomposed into two-particle terms, whose color and momentum structures resemble 
that of one-loop diagrams. At two-loop order, these cancellations were found by explicit 
calculation in [30, 31]. More recently, the analysis was extended to the subclass of three- 
loop graphs containing fermion loops [32]. In Section 6.2 we will present a simple symmetry 
argument explaining these results. 

To derive the perturbative expansion of the Z^-factor from the formal solution (6) we use 
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the generalized expression 



dag 
din// 



= P{as, e) = P{as) - 2e as 



(8) 



for the /3-function in d = 4 — 2e dimensions, where ag = Q!s(Ai) is the renormalized coupling 
constant. The simple form of (7) implies that the matrix structure of the anomalous dimension 
is the same at all scales, i.e., [r({p}, //i), r({p}, //2)] = 0. The path-ordering symbol can 
thus be dropped in (6), and we can directly obtain an expression for the logarithm of the 
renormalization factor. Writing r({j9}, fi, Q;s(/i)) instead of r({p}, //) to distinguish the explicit 
scale dependence from the implicit one induced via the running coupling, we obtain 



as r- a 

, ^ , , , , f da 1 -r,/ r 1 \ f da' r'fo;') 



(9) 



where ag = Q;s(/x), and we have defined 

d 



r'iag) 



dlafj, 



T{{p},i^,as) = -7cusp(as) ^1 



(10) 



Note that this is a momentum-independent function, which is diagonal in color space. We 
have used that, when acting on color-singlet states, the unweighted sum over color generators 
can be simplified, because relation (3) implies that 



(11) 



This relation can be used in our case, because the scattering amplitudes are color conserving. 
Note that a difi^erent but equivalent form of relation (9) has been given in [3] . 

It is understood that the result (9) must be expanded in powers of ag with e treated as a 
fixed 0{a^) quantity. Up to three-loop order this yields 



InZ 



r' 



+ 



An \ 4e2 2e 



16e2 



£1 

4e 



+ 



An) 



72e* 72e3 36e^ 6e 



(12) 

+ 0{a% 



where we have expanded the anomalous dimensions and /3-function as 

r^Er.(S)"". r-f;r^(S)"". 



n=0 



n=0 



n=0 



(13) 



Exponentiating the result (12) and taking into account that the different expansion coefficients 

r„ commute, it is straightforward to derive an explicit expression for Z. For the convenience 
of the reader, we present the result along with the relevant expansion coefficients of the 
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anomalous dimensions in Appendix A. Note that the highest pole in the 0(q;") term of In Z 
is l/e"+^, instead of for the Z-factor itself. The exponentiation of the higher pole terms 
was observed previously in [33]. 

The IR singularities of two-loop scattering amplitudes were first predicted by Catani a 
decade ago [20] . The one- and two- loop coefficients of our Z-matrix are closely related to his 
subtraction operators /^^^ and I^^^ . Catani's formula states that the product 



l-Mn(e, W)) 



(14) 



is free of IR poles through OioPg). The subtraction operators I^^^ie) = /("^(e, {p}, //) are 
defined as 



jW(e) = 



r(l-e) 4-^62 2eT^Jj^^ 2 



(15) 



= (^ + I) ^'"(2.) - ii<'>(.) (/'■'(.) + f ) + HSkie) . 

The conditions linking these objects to ours are 

2J(^) = Zi + finite , 4/(2) = Z2 - 21^^^ Zi + finite , (16) 

where Z„ denotes the coefficient of (as/47r)" in the Z-factor. The first relation is indeed 
satisfied. The second one can be used to derive an explicit expression for the quantity H^'^ 
encoding the genuine two-loop coefficient of the 1/e pole in (15), which was not obtained in 
[20]. We find 



16e ^ 



'Sjk 



In 



;i7) 



Sjk 



'Ski Sij 



which apart from the last term is diagonal in color space and universal in the sense that it 
is a sum over contributions from each individual parton. Note that only the first term in 
this result is of a form suggested by (9). The remaining terms in the first line arise because 
the two-loop corrections involving the cusp anomalous dimension or the /9-function are not 
implemented in an optimal way in (15). More importantly, the term in the second line of (17) 
arises only because the operator I^^^ in [20] is not defined in a minimal subtraction scheme, 
but also includes 0(e") terms with n>0. As a result, the antisymmetric terms in the product 
I^^^Zi in the second relation in (16) contain the structure 



In^ 



-Ski 



Tj, Tk ■ Ti] , 



(18) 
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which after some algebraic simphfications reduces to the expression shown in the second hne 
of (17). Our result for H^^ agrees with the findings of [31] and confirms a conjecture for 
the form of -ff^s ^'^^ ^ general n-parton amplitude made in [34]. Note that the last term in 
(17) is only non-zero for four or more partons. Due to color conservation the three-parton 
case involves only two independent generators Ti and T2, which is not sufficient to obtain a 
completely antisymmetric structure to contract with /"^'^ in (17) or, equivalently, to get an 
non-zero commutator term in (18). In Section 6.2 we will show that the non-trivial color and 
momentum structure in the second line (17) is incompatible with constraints derived from soft- 
collinear factorization, and thus it cannot appear at any order in the anomalous-dimension 
matrix, from which the .Z-factor is derived. 

Our expressions (12) and (17) reproduce all known results for the two-loop 1/e" poles of 
on-shell scattering amplitudes in massless QCD. In addition to the on-shcll quark and gluon 
form factors, these include e^e~ qqg [35] as well as all four-point amplitudes of quarks and 
gluons [36-42]. At the three-loop level, only the IR divergences of the quark and gluon form 
factors are known for the QCD case [26-28]. For A/" = 4 SYM in the planar limit, on the other 
hand, the four-point functions are known up to four-loop order [43, 43], and they also agree 
with our result. 

An interesting alternative approach to the problem of IR singularities of on-shell amplitudes 
was developed in [33], where the authors exploited the factorization properties of scattering 
amplitudes [45-47] along with IR evolution equations familiar from the analysis of the Su- 
dakov form factor [48]. They recovered Catani's result (15) at two-loop order and related the 
coefficient of the unspecified 1/e pole term to a soft anomalous-dimension matrix, which was 
unknown at the time. They also explained how their method could be extended beyond two- 
loop order. The two-loop soft anomalous-dimension matrix was later calculated in [30, 31]. 
In very recent work, Gardi and Magnea have pushed this approach further and derived a set 
of constraint relations for the soft-anomalous dimension matrix, which hold to all orders in 
perturbation theory [29]. We will comment later on the relations between our analysis and 
their work. 

3 Anomalous dimensions of n-jet SCET operators 
3.1 Basic elements of SCET 

To analyze a hard-scattering process involving energetic particles propagating along the direc- 
tions of unit three- vectors hi in SCET, we introduce two light- like reference vectors nj = (1, ri^) 
and fii = (1, —hi) for each direction, so that rii ■ fii = 2. The effective theory then contains 
a set of coUinear quark and gluon fields for each direction of large momentum flow. These 
describe partons with large energies Ei ~ V§ associated with a jet of small invariant mass M. 
The small ratio of these scales, A = M'^/s, serves as the expansion parameter of the effective 
theory. The components of the momenta pc of the coUinear quark and gluon fields Xi ■^i 
associated with the i-th jet direction scale as 

i-coUinear: {fii ■ Pc, rii ■ Pc, p^x) ~ (1; "V^) , (19) 
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such that ~ A ~ M^. The ± components are defined to be perpendicular to both rii and 
Hi. Via the equation of motion, the scahng of the momenta also implies a scaling for the spin 
components of the fields. In the case of collinear fermions, it implies that two of the four 
components of the Dirac spinor field are power suppressed. These can be integrated out, after 
which the field fulfills the condition ^^Xi — 0- Fo^" the collinear gluon field, it implies that 
leading-power operators only depend on Afj_. 

In the absence of soft interactions, each collinear sector of the theory is equivalent to the 
original QCD Lagrangian [7]. This is not surprising, since we can imagine going into the rest 
frame of any given jet, and the interactions among the partons of the jet will then be the 
same as in ordinary QCD. The particular scaling of the fields (19) is dictated by external 
kinematics, or more concretely by the source terms which generate them. The purely collinear 
SCET Lagrangian is thus simply given by n copies of the ordinary QCD Lagrangian, and the 
effective-theory fields Xi '^i± ^-^e related to the usual quark and gluon fields via 

X,{x) = Wl{x)^Aix), Alix) = Wl{x)[tD^^W,{x)]. (20) 
The i-coUinear Wilson lines 

Wi{x) — P exp J dsfii ■ Ai{x + sfii)^ (21) 

ensure that these fields arc invariant under collinear gauge transformations in each sector [4, 5] . 
The symbol P indicates path ordering, and the conjugate Wilson line W- is defined with the 
opposite ordering prescription. 

In addition to coUinear fields for each direction, the effective theory contains a single set 
of soft fields, which interact with all types of collinear fields. All components of the momenta 
Ps carried by these fields scale as 

soft: ~ A . (22) 

This scaling is such that one can associate a soft parton with any of the n jets without 
parametrically changing the invariant mass of the jet, because (pc + Ps)^ ~ A. The soft 
fields can thus mediate low-energy interactions between different collinear fields. However, 
at leading power this interaction is very simple: it can be obtained from the substitution 
Af(x)^Af(x) + fni-A(x_),or 

A^{x) ^ A^ix) + ^ W^{x)gn, ■ A{x_) W,{x) , (23) 

in each of the collinear Lagrangians, where a;_ = (n^ ■ x) ni/2. Only the rii ■ As component of 
the soft gluon field enters in this relation, since all other components are power suppressed 
compared to the collinear fields. The peculiar x-dependence of the gluon field is a consequence 
of the multipole expansion [7, 49], which implies that in interactions of collinear and soft fields 
one should perform a derivative expansion of the form 

[Mx)f <i>six) = [<j),{x)f [(j)s{x_) +x^-d^ (Psix) U=,_ + ...]. (24) 
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The first-derivative term is suppressed by -\/A, because ~ 1/ -\/A, while d± acting on a soft 
field counts as 0{\). 

The substitution (23) gives rise to an eikonal interaction of soft gluons with coUinear 
fermion fields, 

J^ci+s = Xi{^) y • A{x-) Xi{x) ■ (25) 

This interaction can be represented in terms of soft Wilson lines. Redefining the quark and 
gluon fields as 

Xi{x) = Si{x_)xf\x) , 

U^)^xP{x)S}{x^), (26) 



AtA^) = S,{x.)A';Ax)Slix.), 



where 



Si{x) = P exp (^g J dtUi- A'^^ix + tni)t''^ , (27) 

eliminates the interaction Ca+s (including the pure-gluon terms). After this decoupling trans- 
formation [5], soft interactions manifest themselves as Wilson lines in operators built from 
collinear fields. The soft gluons do not couple to the spin of the collinear particles, and for 
the discussion that follows the spin degrees of freedom will be irrelevant. 

As written above the soft Wilson lines Si and Sj are color matrices defined in the fun- 
damental representation of the gauge group. The transformations (26) take on a universal 
form if wc define a soft Wilson line Si in analogy with (27), but with t° replaced by the color 
generator in the appropriate representation for the i-th parton. Representing a generic 
collinear field as {4>i)2l{x) with color index aj and Dirac/Lorentz index ai, the soft interactions 
can then be decoupled from this field by the redefinition 

{<p^::i^)-[s^{x-)UA<p^f^::i^)■ (28) 

Note that even anti-quarks transform according to this rule: in this case — — (t")^, which 
translates into the anti-path ordering in (26). 

Hard interactions among the different jets are integrated out in the effective theory and 
absorbed into the Wilson coefficients of operators composed out of products of collinear and 
soft fields. Since additional soft fields in the SCET operators would lead to power suppression, 
the leading n-jet operators are built from n collinear fields, one for each direction of large 
energy fiow [50, 51]. The most general such operator with given particle content appears in 
the effective Hamiltonian 

nf = Jdh... dtnCZ-:X{h, ...,tn,l^) {(j)iral{x + hu,) . . . {(j)n)Zi^ + ^n^n) • (29) 

Our notation is somewhat unusual, because the Wilson coefficients of these operators carry 
spin and color indices. Usually both operators and Wilson coefficients are chosen to be color- 
neutral Lorentz scalars. However, writing the operator in this form makes the connection to 
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the color-space notation we use for the scattering amphtudes most transparent. Adopting this 
notation, the effective Hamiltonian for an n-jet process reads 

nf = Jdt,... dt^ (o„({t},/.)|c„(U},/x)) , (30) 

where ji is the scale at which the SCET operator is renormalized. An n-parton scattering 
amplitude is obtained by taking an on-shell partonic matrix element of this operator. In 
this step effective-theory loop integrals vanish in dimensional regularization, because they are 
scaleless. The on-shell matrix elements are therefore given by their tree-level values, consisting 
of products of on-shell spinors and polarization vectors defined through the relations 

(0| (Xj)a(^j%) \Pi] ai, Si) = 5ij 5a,a e-'^'^^'P' Ua(Pi, Si) , 

(31) 

{^\i^j±)t{tjnj) \Pi] ai, Si) = 6ij 5a,a e e^(pj, Sj) . 

Loop corrections to the scattering amplitude are encoded in the Wilson coefficients C„({t}, n). 
The integrations over tj in (30) produce the Fourier transforms C„({p}, ^) of these coefficients, 
which after contraction with the spinors and polarization vectors arising when taking the tree- 
level matrix elements are in one-to-one correspondence with the scattering amplitudes [3], as 
shown in (2). 



3.2 Soft-coUinear factorization and decoupling 

To obtain the general form of the anomalous-dimension matrix V defined in (4), we now derive 
a factorization theorem for the matrix elements of SCET operators. The first factorization step 
has already been achieved in (30), which separates hard from soft and collinear fluctuations. 
In a second step, we separate the collinear and soft degrees of freedom using the decoupling 
transformation (26), which eliminates the leading-power interactions among soft and collinear 
fields. Since collinear fields from different sectors do not interact directly, this completely 
factorizes a matrix element into a soft part 5, given by a product of Wilson lines along the 
directions n^, and a product of collinear matrix elements Ji for each direction. 

RG invariance implies that the right-hand side of (30) must be independent of the rcnormal- 
ization scale. Denoting by F/j = F the anomalous-dimension matrix of the hard contributions 
contained in the Wilson coefficient functions C„ and by Fc+s the anomalous dimension associ- 
ated with the collinear and soft contributions contained in the matrix elements of the SCET 
operators,^ it follows that = F^+s. The decoupling transformation, which removes the 
interactions of collinear fields with soft gluons and absorbs them into Wilson lines [5], allows 
us to further decompose Fg+s = Fg + F^. There are no mixed soft-coUinear contributions. The 
collinear piece Fc = F* is a sum over color-singlet single-particle contributions, because the 
fields belonging to different collinear sectors of SCET do not interact with one another. Hence, 
contributions to the anomalous dimension involving correlations between different partons only 
reside in the soft and hard contributions, F^ and T^, and they coincide. 

^Following common practice, we define the anomalous dimensions of operators with the opposite sign 
compared to those for Wilson coefficients. 
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Figure 1: Graphical representation of the soft operator S{{n}, fi) corresponding to an n-parton 
scattering amphtude. The n hght-like Wilson lines start at the origin and run to infinity. The 
dots represent open color indices. 

After the decoupling transformation the soft matrix element is a vacuum expectation value 
of n light-like Wilson lines, one for each external parton in the associated color representation: 

5({z^},/^) = (0|5i(0)...5„(0)|0). (32) 

As illustrated in Figure 1, this object is an operator in color space, with each Si factor 
operating on the color indices of the i-th parton. Its renormalization properties are strongly 
constrained by the simplicity of soft gluon interactions, which only probe the direction of the 
Wilson lines and their color charge. When the color indices are contracted in color-singlet 
combinations, then «S({n},//) turns into products of closed Wilson loops, which touch or 
intersect each other at the origin. The renormalization properties of such Wilson loops have 
been studied extensively in the literature, see e.g. [8, 9, 23, 52-55] and references therein. 
We will use several results obtained in these studies and generalize them to the case of the 
Wilson-line operator in (32). We will also indicate where known properties of Wilson loops 
correspond to certain features of the effective theory and vice versa. 

For on-shell amplitudes, the loop integrals in the effective theory have both IR and UV 
divergences and vanish in dimensional regularization. This makes the correspondence between 
the Wilson coefficients in (30) and the amplitudes manifest. However, because of the cancel- 
lations between UV and IR poles, we cannot use on-shell amplitudes to obtain the anomalous 
dimensions of the SCET operators. To separate out the UV divergences we need to consider 
IR-finite quantities. The simplest possibility is to consider slightly off-shell n-parton ampu- 
tated Green's functions Gn{{p})- However, in this case we encounter a subtlety. While the 
off-shell Green's function in QCD and SCET are IR finite, this is no longer the case after 
the field redefinition (26). Field redefinitions leave "physical" quantities such as on-shell ma- 
trix elements unchanged, but they can change the off-shell behavior of fields. To calculate the 
anomalous dimensions perturbatively from off-shell Green's functions, one should use the orig- 
inal, non-decoupled fields.^ For the case of the quark form factor, the corresponding one-loop 
calculation in the effective theory was performed in [56]. Generalizing this result to n-point 

^Alternatively, one can perform the calculations using a different IR regulator, e.g. by considering finite- 
length Wilson lines with nf ^ [25]. 
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functions, we find for the UV poles of the jet and soft functions (normahzed to 1 at tree level) 



a. 



47r 



47r 



p 

+ - In- 



2e 



(33) 



2(-pD(-P?) 



+ C»(e°) . 



2 2 , 
^ + - In 



(34) 



The jet functions are color- diagonal. The reference vectors Ui are defined such that up to 
power corrections pi = Eifii, which implies that at leading power | aij rii ■ rij fii ■ pi fij ■ pj — 
'^'^ijPi ' Pj = Sij. The one-loop divergences of the complete effective-theory n-particle matrix 
elements are thus given by 

5(M,/.) n ^^^pit') = 1 - s f 5] ^ 

The UV divergences of the SCET operator matrix elements arc equal and opposite to those of 
the corresponding Wilson coefficient. They are thus given by minus the one-loop coefficient of 
In Z in (12). The one-loop coefficients of the jet-function anomalous dimensions are obtained 
from (33) as 7q — —3Cf and 7q = —/3o- Interestingly, the dependence on the off-shellness 
cancels in the divergent part of the full matrix elements, which only depend on the hard scales 
Sij. This cancellation has to occur, because the anomalous dimensions of the operators cannot 
depend on low-energy scales such as p^. Otherwise the scale dependence of the operators could 
not possibly cancel against the scale dependence of the hard Wilson coefficients, which are 
insensitive to IR scales. As we will explain below, the factorization property F/j = -|- 
places a strong constraint on the form of the anomalous dimension. 



3.3 Color-symmetrized soft gluon attachments 

Of special importance to our analysis in Section 4 will be the fact that the interactions of 
soft gluons with collinear partons take on a particularly simple form after a symmetrization 
of color matrices has been performed. To derive the form of this relation, consider the sums 
over all possible attachments of one or two soft gluons to a set of I Wilson lines Si{0) with 
tangent vectors rij, such as those contained in the soft operator in (32). Using the Feynman 
rules of SCET, we find the simple expressions shown in Figure 2, where the dots in the second 
relation represent an anti-symmetric color structure containing the commutator [T^", J]^] for 
a single parton index. Such a structure can be reduced to a single color generator using the 
Lie algebra of the gauge group. The important point to note about this result is that for the 
symmetric color structure the momentum dependence is the same irrespective of whether the 
soft gluons attach to the same or to different Wilson lines. In the two-gluon case, this happens 
because the two possible attachments to a single Wilson line yield 

rpb rpa rpa rpb 

' i \ i i (•35) 

Hi ■ (ki + /C2) Hi ■k2 Ui - (ki + k2) rii-ki' 
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Figure 2: Feynman rules for the sums over all possible attachments of one (top) or two 
(bottom) soft gluons to a set of Wilson hnes. The indices i,j run from 1 to I. 



which after symmetrization gives rise to the structure shown in Figure 2. 

The Feynman rules given in the figure generalize in an obvious way to the case of more 
than two soft gluons. The reason is that the symmetrization in color generators eliminates the 
need for the path-ordering symbol in the definition (27) of the soft Wilson lines. Introducing 
the notation 

(rpai rpa„\ — ^ \ ^ rp"-a(l} rp"-a(n) ( 'iP.\ 



ni 



for the symmetric product of n color matrices, where Sn is the set of permutations of n objects, 
we obtain for the general case of m soft gluons attached to a set of Wilson lines the Feynman 
rule 



E 



n^} . . . 7if " 



gm ■■■ ir. +.... (37) 

7ljj ■ Ki . . . 7lj^ • Km 



Since color generators acting of different partons commute, the symmetric product in (36) 
symmetrizes the color matrices acting on each individual parton. For instance, we have for 
different indices i, j, k 

(T-T^T^T^T^)^ = {nn)^T^ ■ (38) 



4 Constraints from soft-collinear factorization 

We will now summarize the arguments that have led us to propose the conjecture (7). Ap- 
proaching the problem from the perspective of effective field theory adds one crucial element 
to the discussion. This is the realization that the IR-divcrgent terms in the scattering am- 
plitudes must obey some restrictive constraints, which do not apply to the finite remainders. 
The scattering amplitudes are in general complicated functions of the kinematical invariants 
Sij as well as of the color, spin, and polarization quantum numbers of the external partons. 
The different Sij variables are assumed to be hard scales of the same order of magnitude, so 
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Figure 3: Color-singlet contractions of four Wilson lines in the fundamental representation. 
The resulting Wilson-loop operators mix under renormalization. 

that the ratio of any two such variables is an 0{1) quantity. In principle, arbitrary functions 
of combinations of such ratios can arise in the expressions for the scattering amplitudes. The 
situation is, however, very different for the IR-singular terms in the amplitudes. RG invariance 
of the effective theory requires that the anomalous dimensions of the hard matching coeffi- 
cients /Li)), which according to (2) correspond to the on-shell scattering amplitudes, 
must be decomposable into sums of collinear and soft contributions. This requires a rewriting 
of the hard momentum variables Sij in terms of soft and collinear variables. The very fact that 
such a rewriting must exist restricts the functional dependence of the anomalous dimension 
on the Sij variables to be either single logarithmic or in the form of certain conformal cross 
ratios, which will be defined below. Moreover, the structure of the effective theory enforces 
that terms depending on the collinear variables cannot lead to correlations between different 
partons and must be diagonal in color space. Correlations can only arise through soft gluon 
exchange. The universal structure of these interactions implies that any dependence on the 
identity of the external partons can only arise via their momenta and color charges, but not 
through spin information. We will also discuss constraints on the color structure of the soft 
anomalous-dimension matrix implied by the non-abelian exponentiation theorem and other 
considerations. 

4.1 Renormalization of Wilson loops 

A well-known property of Wilson loops is that they require UV subtractions beyond the renor- 
malization of the coupling constant in cases where the integration path is not smooth, but 
contains one or more singular points [52-54]. These divergences can be removed multiplica- 
tively. The simplest case is that of a Wilson loop with a single cusp, i.e., a point where the 
tangent vector changes its direction abruptly. If the cusp is formed by two time-like segments 
with tangent vectors rii and n2, then these UV divergences are removed by a factor Z{/3i2), 
which is a function of the hyperbolic cusp angle f3i2 defined by 

cosh/3i2 = ^^^, (39) 

where for simplicity we have assumed that rii points into the cusp and n2 points out of it. If 
the Wilson loop has more than one cusp, then each of them contributes an analogous Z-factor. 

A more complicated situation arises if, as in our case, different Wilson lines cross each 
other at a point. Then Wilson loops tracing out the same space-time curves except for the 
cross point mix under renormalization. An example are the two Wilson-loop operators shown 
in Figure 3, which illustrates this fact for the case of a four-jet operator corresponding to 
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Figure 4: Diagrammatic form of the Lie commutator relation. Gluons are drawn as wavy lines 
in order to distinguish color- weight graphs from Feynman diagrams. 

qq — > qq scattering. The renormalization factor Z({/3}) is then a matrix on the space of such 
Wilson loops, which depends on the set of all hyperbolic angles formed by the tangent vectors 
at the cross point [54]. Generalizing these results to our case, where the Wilson-line operators 
are matrices in color space as shown in (32), the renormalization factor must be promoted to 
a soft matrix Zg acting on the product space of the color representations of the n partons. 

RG invariance implies that the renormalization factor can be constructed in the usual 
way from a soft anomalous-dimension matrix F,. For the case of a single cusp, the two-loop 
expression for the anomalous dimension was first obtained in [23, 24]. 

4.2 Non-abelian exponentiation theorem 

The structure of the soft anomalous-dimension matrix is restricted by the non-abelian ex- 
ponentiation theorem [8, 9], which implies that purely virtual amplitudes in the eikonal ap- 
proximation (i.e., with only soft gluon interactions taken into account) can be written as 
exponentials of simpler quantities, which only receive contributions from Feynman graphs 
whose color weights are "maximally non-abelian" (or "color-connected"). Applied to our case, 
it follows that the logarithm of the soft Z^-factor, and with it the soft anomalous-dimension 
matrix F^, only receives such contributions. 

We follow the diagrammatic approach to the non-abelian exponentiation theorem devel- 
oped in [9], since it is more explicit and intuitive than the iterative construction presented in 
[8] . To each Feynman diagram we assign a color- weight diagram, in which vertices are replaced 
by color matrices {t"')ij or structure constants —if"-^'^ (or, more generally, by generators T"" 
in the appropriate representation of the gauge group), and propagators by ^'-^ for quarks and 
^ab gluons or ghosts. Color diagrams may be related to one another by use of the Lie 
algebra relation [T", T^] = if"-^'^ T^, as illustrated in Figure 4. In the adjoint representation 
this is called the Jacobi identity. One defines a connected web as a connected set of gluon lines 
(possibly with internal fermion loops), not counting crossed lines as being connected. It has 
been shown in [9] that, using the Lie commutator relation, any color- weight diagram can be 
written as a sum over products of connected webs. The non-abelian exponentiation theorem 
implies that only single connected webs contribute to the color weights in the exponent. For 
the case of two Wilson lines this is illustrated with an example in Figure 5. In this special 
case, the above definitions imply that the single connected webs contain those diagrams that 
are two-particle (rainbow) irreducible diagrams with respect to the Wilson lines [23]. 

Note that in our case the gluons can be attached to more than two Wilson lines, provided 
there are more than two external partons. The notion of crossed lines is meaningful only for 
diagrams such as Figure 5, in which only two Wilson lines are connected. The reduction to 
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Figure 5: Decomposition of a web into a sum of products of connected webs. The non-abelian 
exponentiation theorem states that only the single connected web shown in the first graph on 
the right contributes to the color weights in the exponent of the amplitude. 

connected webs needs to be generalized when more than two lines are connected by gluons in 
a diagram. In general, if one has a diagram with n gluons connected to a given Wilson line i, 
one first symmetrizes the product of color matrices using the identity 

T^K.. T^^ = {T^^ . . . i;"")^ + "commutator terms" . (40) 

The "commutator terms" involve one or more commutators of color generators and can thus 
be reduced to products of (n — 1) or fewer generators multiplied by structure constants. Re- 
peatedly applying the identity (40), one can write any color structure as symmetric products 
of generators multiplied by structure constants. By performing the symmetrization on all legs, 
one splits each diagram into a number of webs, and the statement of non-abelian exponentia- 
tion is that only single connected webs contribute in the exponent. Since all attachments are 
fully symmetrized, the distinction between crossed and uncrossed lines is no longer relevant. 

The fact that only single connected webs contribute to the logarithm of the .Z^-factor 
(and hence to the anomalous dimension), while products of webs contribute to the Z^-factor 
itself, is in analogy to the usual structure of UV divergences in quantum field theory, as 
described by Zimmermann's forest formula. Formal arguments explaining the systematics of 
UV divergences for arbitrary Wilson loops can be found in [54] . 



4.3 Light-like Wilson lines 

For large values of the cusp angle /3i2 in (39), the anomalous dimension r(/3i2) associated with 
a cusp (or cross point) grows linearly with (5i2 [23], which in this case is approximately equal to 
ln(2ni-n2/ \/ n\ n\). In the limit where one or both segments forming the cusp approach a light- 
like direction, the cusp angle diverges (/3i2 oo). In dimensional regularization this divergence 
gives rise to a single logarithm of the renormalization scale in the anomalous dimension. If 
both segments he on the hght-cone, then [25] 

r(A2) VI (a,) In ^ + . . . , (41) 
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where we refer to r*^gp(Q;s) as the cusp anomalous dimension in the color representation of 
parton i. Its two-loop expression was obtained long ago in [57, 58] and [23, 24], while the 
three-loop result was derived in [13]. The above equation is formal and meant to show the 
dependence on the renormalization scale only. We will explain later how a soft scale with 
the proper dimensions appears in the argument of the logarithm. 

In conventional applications of the RG, large (single) logarithms of scale ratios entering 
perturbative results for multi-scale problems can be resummed with the help of anomalous 
dimensions that are functions of the coupling constant, much like the /3-function. This resums 
terms of the form (a^L)" in the perturbative series, where L is the logarithm of the relevant 
scale ratio. However, the presence of overlapping soft and coUinear singularities in on-shell 
scattering amplitudes of massless partons generates Sudakov double logarithms of the form 
(ttsL^)" in perturbation theory. They can be resummed with the help of anomalous dimen- 
sions which themselves contain a single logarithm L of the large scale ratio. The logarithmic 
dependence of the anomalous dimension in (41) is an essential feature in this context. 

SCET is the appropriate effective field theory to formalize apphcations of the RG to ob- 
servables sensitive to Sudakov double logarithms. It has been used in the past to derive the 
exact form of the anomalous dimensions of the two-jet operators for quarks and gluons, which 
were found to be of the form [59-62] 

r^-jet = -K^^Ms) 1^ + 2f («.) ■ (42) 

Here i — q for quarks and i = g for gluons, and s = Si2 is the invariant momentum transfer. 
From the perspective of the effective theory, the appearance of the cusp logarithm is, at first 
sight, perplexing. How can the anomalous dimension of an operator defined in an effective 
theory, in which the hard scale s has been integrated out, remember the value of that scale? 
The answer to this puzzle was given in [56], where it was explained that the hard scale is 
imprinted in the effective theory via the large rapidity that separates the rest frames of soft 
and coUinear hadrons in a given physical process. This leads to a characteristic entanglement 
of the hard, coUinear, and soft mass scales, the latter two of which are known to the effective 
theory. This correlation is such that /i^ ~ l^^hfJ^s, and the hard logarithm in relation (42) can 
thus be rearranged in the symbolic form 

ln4 = 21n^-ln^. (43) 

Obviously, such a rearrangement of a hard contribution as a sum of coUinear and soft con- 
tributions is only possible for functions containing either constants or single logarithms of 
scale ratios, a point that was also emphasized in [62]. We observe an interesting connection 
between RG invariancc in SCET and a property of soft Wilson loops: RG invariance requires 
a single-logarithmic dependence of the anomalous dimension on the hard scale, because only 
then can this dependence be decomposed into dependences on coUinear and soft scales. Re- 
lation (43) then implies single-logarithmic dependence on the soft scale, in accordance with 
known renormalization properties of Wilson loops mentioned above, see (41). 

That a decomposition of the form (43) is indeed at work in the effective theory was demon- 
strated in [56] using the method of regions, by analyzing the coUinear and soft contributions 
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to the anomalous dimension (42) separately. As we have explained in Section 3.2, in doing this 
it is necessary to either consider physical quantities, which are IR finite, or to introduce an 
IR regulator scale in order to define the collinear and soft scales. Using a regulator introduces 
some arbitrariness and scheme dependence into the calculation of the individual contributions, 
which however vanishes in their sum (see also [29, 31]). For concreteness, we introduce a small 
off-shellness (— pf) > for the external partons, taking the hmit pf — > wherever possible. 
The decomposition of a generic hard logarithm then reads 

In^^ln ^ = lni^ + ln^-ln-'^-^^;-^^/. (44) 



-s 



'i'JijPi-Pj -pi -p] (-pf)(-p?) 



13 ^^i]fi fj fi fj 

This is precisely the structure of collinear and soft logarithms found in [56]. Measuring all 
scales in units of the hard scale, we have the power counting pi ■ ~ 1 for the hard scales, 

~ p| ~ A for the collinear scales, and pfp^/pi ■ Pj ~ for the soft scales, in accordance 
with the general discussion in Section 3.1. 

In our discussion below we will assume that such a regularization is employed. We then 
introduce the notations^ 

{-Pi){-Pj) -Pi 

for the logarithms of the soft and collinear scales, respectively. The definition of Pij generalizes 
that of the cusp angle in (39) to the case of light-like Wilson lines. The role of the soft scale 
in (41) is played by = (— p|)(— p^)/(— 2cryPj -Pj). Relation (44) can now be rewritten as 

i3ij = Li + Lj-]n-^. (46) 

Sij 

4.4 General structure of the soft anomalous-dimension matrix 

We are now ready to analyze the structure of the anomalous-dimension matrix of n-jct SCET 
operators. According to the discussion in Section 3.2, we can write the decomposition into 
soft and collinear pieces as 

r({p}. Pi) = r,({^}, + KiLi: Pi) , (47) 

i 

where the collinear terms are diagonal in color space. The total anomalous dimension depends 
on the n(n — 1)/2 kinematical variables Sij, while its soft counterpart depends on the n(n— 1)/2 
cusp angles f3ij, as indicated. The collinear pieces are single-particle terms, each depending on 
a single collinear scale Li . The general form of the collinear part of the anomalous dimension 
is known to be of the form [56] 

r:.(L,) = -n„,p(«.)L, + 7>,). (48) 



^At leading power in the effective theory, the product 2pi ■ pj in the argument of the first logarithm is 
replaced by ^rti- nj fii ■ pi fij ■ pj, see (33). 



20 



We now substitute for the cusp angles entering the soft anomalous-dimension matrix the 
expression on the right-hand side of (46). This yields ^s{{§}, {L}, fi) as a function of the 
variables Sij and Lj. The dependence on the collinear scales must cancel when we combine 
the soft and collinear contributions to the total anomalous-dimension matrix. We thus obtain 
the relation 

'-^^f^-TU^.). (49) 

where the expression on the right-hand side is a unit matrix in color space. This relation 
provides an important constraint on the momentum and color structures that can appear in the 
soft anomalous-dimension matrix. A corresponding relation has been derived independently 
in [29]. 

Because the kinematical invariants Sij can be assumed to be linearly independent, relation 
(49) implies that Ts depends only linearly on the cusp angles (3ij, see (46). The only exception 
would be a more complicated dependence on combinations of cusp angles, in which the collinear 
logarithms cancel. The simplest such combination is 

= Pij + Pki - Pik - Pji = In [Z'^'WZ"''] , (50) 

[ Sik)[ Sji) 

which coincides with the logarithm of the conformal cross ratio pijki defined in [29]. For 
simplicity, we will use the term "conformal cross ratio" in the following also when referring to 
Pijki- This quantity obeys the symmetry properties 

Pijkl — Pjilk — —Pikjl — —Pljki — Pklij ■ (51) 

It is easy to show that any combination of cusp angles that is independent of collinear loga- 
rithms can be expressed via such cross ratios. Moreover, given four parton momenta, there 
exist two linearly independent conformal cross ratios, since 

f^ijkl + Awj + f^iljk = , (52) 

and all other index permutations can be obtained using the symmetry properties in (51). 

Our strategy in Section 6 will be to analyze the structure of the soft anomalous-dimension 
matrix first, since it is constrained by the non-abelian exponentiation theorem and the con- 
straint (49). The universality of soft gluon interactions implies that the soft contributions only 
probe the momentum directions and color charges of the external partons, but not their po- 
larization states. Dependence on the parton identities thus only enters via the cusp variables 
f3ij and non-trivial color-conserving structures built out of Tj generators. If our conjecture (7) 
is correct, then (48) implies that the soft anomalous- dimension matrix should be given by 

where 

^\as) = I'Ms) + il{as) . (54) 
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Using relation (11) we may indeed confirm that 
dT 

-qJ^ ^ - ^^i-^j 7cusp(Qis) = Ci 7cusp(Qi.) = KnspM , (55) 

in accordance with the constraint (49). Note that this result implies Casimir-scaling for the 
cusp anomalous dimension, since T^^^^^ag) /T1^^gp{as) — Ca/Cf- We will come back to the 
significance of this observation in Section 6.4. 



5 Consistency with collinear limits 

Before turning to a diagrammatic study of the anomalous-dimension matrix we discuss one 
more non-trivial constraint it must obey, which derives from the known behavior of scattering 
amplitudes in the limit where two or more external partons become collinear. 

In the limit where the momenta of two of the external partons become collinear, an n-parton 
scattering amplitude factorizes into the product of an (n — l)-parton scattering amplitude 
times a universal, process-independent splitting amplitude. This was first shown at tree level 
in [63, 64], and extended to one-loop order in [65]. An all-order proof was given in [66]. 
Strictly speaking, the proof was constructed for leading-color amplitudes only, but the crucial 
ingredients are unitarity and analyticity, and it should be possible to extend it to the general 
case. Collinear factorization holds at the level of the leading singular terms. It is often studied 
for color-ordered amplitudes, for which the color information is stripped off. The color-stripped 
splitting amplitudes for the splitting of a parent parton P into collinear partons a and h are 
usually denoted by Split^^(a'^", 6'^*') in the literature, where cXi denote the helicities of the 
partons. These functions have been calculated at tree level (see, e.g., [67]), one- loop order [68], 
and recently even to two loops [69]. In contrast, we will study collinear factorization using the 
color-space formalism, extending the work of [70] beyond the one-loop approximation. In this 
framework, the splitting amplitudes are elements of a splitting matrix Sp({pa,pb}), which acts 
in the space of color and spin configurations of (n — l)-parton scattering amplitudes. As is the 
case for the scattering amplitudes, the divergence structure of Sp({pa,p5}) is independent of 
the spin configuration of the involved partons, and we therefore suppress spin indices in the 
following. For Catani's formula (15), the consistency with collinear limits was shown in [34]. 

Consider, for concreteness, the limit where the partons 1 and 2 become collinear and merge 
into an unresolved parton P. We assign momenta p\ — zP and p2 = {1 — z)P and consider 
the collinear limit ^ 0. In this limit the scattering amplitude factorizes in the form 

\Mn{{Pl:P2:P3: ■ ■ ■ :Pn})) = Sp({pi,P2}) \Mn-l{{P:P3: ■ ■ ■ :Pn})) + ■■■ ■ (56) 

The matrix of splitting amplitudes encodes the singular behavior of the amplitude \Ain) 
as Pi||p25 and the factorization holds up to terms that are regular in the collinear limit. 
Analogous relations describe the behavior in limits where more than two partons become 
collinear. However, it is sufficient for our purposes to focus on the simplest case. 

The factorization formula (56) holds both for the dimensionally regularized scattering 
amplitudes |A^„(e, {p})) as well as for the minimally subtracted amplitudes |A1n({p}, A*)) in 
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(1). Since we know that the divergences of the amphtude can be absorbed into a Z-factor, 
equation (56) imphes a constraint on the divergences of the sphtting amphtudes. It can be 
written as 

hniZ-^(e, {pi, . . . ,pn}, n) Sp(e, {pi,P2}) Z{e, {P,p3 . . . = Sp({pi,p2}, /^) , (57) 

where the matrix of renormahzed sphtting amphtudes on the right-hand side is finite for e ^ 0. 
From (5) it then foUows that this quantity obeys the RG equation 

Sp({pi,p2},/^) =r({pi,...,p„},/x)Sp({pi,p2},/^) 



- Sp({pi, P2}, /^) r({P, P3...,pn},ll). 

Analogous equations hold for the higher splitting amplitudes Sp({pi, . . . A*); which de- 
scribe the limits where more than two partons become collinear. To bring the RG equation 
into a more useful form, we note that charge conservation implies 

(Ti + T2) Sp({pi,p2}, 1^) = Sp({pi,P2}, 1^) Tp , (59) 

where Tp is the color generator associated with the parent parton P. Since the matrix of 
splitting amplitudes commutes with the generators of partons not involved in the splitting 
process, we can thus commute the anomalous dimension in the second term to the left to 
obtain ^ 

^|^Sp({pi,p2},/^) = rsp({pi,p2},/^) Sp({pi,p2},/x) , (60) 
where we have defined 

rsp({pi,P2},/^) = r({pi, . . . ,p„},/i) - r{{p,p3 . . . ,p„},/i)|^^^^^_^^^ . (61) 

The fact that the anomalous dimension of the splitting amplitudes must be independent 
of the colors and momenta of the partons not involved in the splitting process, which is a 
consequence of the factorization formula (56), imposes a non-trivial constraint on the form of 
the anomalous-dimension matrix. We will explore its implications in Section 6.6. 

Assuming the form (7) for the anomalous-dimension matrix F, we find that the anomalous 
dimension of the splitting amplitudes has the simple all-order form 



Ti • T2 In + Ti • (Ti + T2) In^ + T2 • (Ti + T2) ln(l - z) 
—S12 



rSp({Pl,P2},/^) = 7cusp 

+ 7^+72-7^, (62) 



where 'y^ is the anomalous dimension associated with the unresolved parton P. Note that the 
momentum-dependent terms in the result are insensitive to the nature of the partons involved 
in the splitting process. The divergent part of the one-loop splitting amplitudes for m partons 
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in the color-space formalism was given in [70]. Expanding the result obtained there for the 
case m = 2, we find 



SPi-ioop(e> {PuP2}) = -r- 



47r 
2 

+ - 
e 



2 2 ^2 
— + - In ■ 



-Sl2 



Ti • (Ti + T2) In ^ + T2 • (Ti + T2) ln(l - z) 

SPtree({Pl>P2}) , 



(63) 



which is in agreement with the result obtained by solving the RG equation (60). Beyond one- 
loop order the splitting amplitudes are given by rather complicated expressions [69] ; however, 
we have checked that their IR divergences can indeed be obtained from the simple anomalous 
dimension in (62), which only contains single logarithms of the momentum fractions z and 
(l — z). The exact form of the anomalous-dimension matrix in (62) is an important by-product 
of our analysis. 



6 Diagrammatic analysis 

We now present a detailed diagrammatic study of the general structure of the soft anomalous- 
dimension matrix up to three-loop order, implementing the constraints that follow from the 
non-abelian exponentiation theorem and RG invariance of the effective theory. At two-loop 
order we will recover the form found in [30, 31] from a simple symmetry argument. In these 
papers only the cusp piece of the soft anomalous-dimension matrix was studied, which is 
legitimate given that the non-logarithmic terms can be shown to be diagonal in color space. 
We find that this property is no longer trivial beyond two-loop order. 

The non-abelian exponentiation theorem restricts the color structures that can potentially 
appear in the soft anomalous-dimension matrix. They are obtained by considering single 
connected webs, whose ends can be attached in arbitrary ways to the n Wilson lines in the 
soft operator in (32). In general, single connected webs at L-loop order carry between 2 and 
(L + 1) color generators T. In Figure 6 we show the webs appearing up to three-loop order. 
The dashed blobs represent self-energy or vertex functions, which have color structure 5"^ and 
_j^jobc fpjjg color structures of the three- and four-gluon vertices can be expressed in terms 
of /"^'^ symbols. 

In our analysis in this section we only use basic properties of the Lie algebra of the gauge 
group, which can be summarized as 

I'/yia ^ jabc rj-ic jabc jabd C j\ 5'^ 

jC (64) 

t-'^ad' {T"" X"^) i J"''^'^ J^^9 j<^9'l ''^ A jabc 

The last relation follows from the Jacobi identity, i.e., the first relation in the adjoint repre- 
sentation. While our explicit analysis refers to SU (N) non-abelian gauge theories, its validity 
extends to other gauge groups as well. 
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Figure 6: One- loop (a), two-loop (b), and three-loop (c) connected webs contributing to the 
soft anomalous- dimension matrix. The dots represent color generators, which appear when 
the gluons are attached to Wilson hnes. In each set, only the first web gives rise to a new 
color structure. 

6.1 One-loop analysis 

In this case the relevant web consists of a single gluon, as shown in Figure 6(a). If it is attached 
to two different Wilson lines i and j, then the resulting color structure is Tj • Tj. In this case 
non-trivial momentum dependence can arise, which can lead to a factor Recall that only 
linear dependence on the cusp angle is allowed. For terms without momentum dependence, 
the sum over parton legs reduces the color structure to a diagonal one, since relation (11) can 
be applied in this case. Likewise, if the ends of the exchanged gluon are attached to a single 
Wilson line i, then the color structure is = Ci. It follows that at one-loop order the soft 
anomalous-dimension matrix is indeed of the form (53). 

6.2 Two- loop analysis 

In this case two webs need to be considered, which are depicted in Figure 6(b). The connected 
web containing the gluon self-energy has the same color structure as a single gluon exchange, 
and hence it does not lead to any new structures in the result (53). The color structure of 
the three-gluon web is proportional to —if'^^'^ times three color generators, one for each leg. 
There are thus three possibilities, which we consider separately. 

If all gluons are attached to a single Wilson line, then the resulting color structure is 



In this case no momentum dependence can arise. If the gluons are attached to two different 
Wilson lines i and j, then the resulting color structure is (recall that generators belonging to 
different partons commute) 




CaCi 



(65) 



2 




(66) 
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In this case momentum dependence can arise, since two partons are involved in the loop 

diagram. It is thus possible to get a factor f3ij or a constant. In any case we obtain the same 
structures as at one-loop order. Finally, if the gluons are attached to three different Wilson 
lines i, j, and k, then the resulting color structure 

is totally anti-symmetric in the parton indices, and it would therefore need to multiply a 
totally anti-symmetric momentum-dependent structure formed out of the three kinematical 
invariants Pij, Pjk, and Pki- However, no such structure exists that would be consistent with 
our constraint (49), since it requires linearity in the cusp angles in cases where less than four 
parton momenta are involved. A nonlinear structure such as 

{Pij — Pjk){Pjk — Pki){Pki — Pij) 

^(l^Zfil + Li- Lk] (in ^ + Lj - L,) (in ^ + Lk - Lj) , ^^^^ 
V -Sjk J V -Ski J V -Sij J 

which is reminiscent of that appearing in (17), cannot be written as a sum of hard and 
coUinear contributions. Recall that this structure arises in Catani's subtraction operator J^^^ 
only because in [20] the subtraction of IR poles is not implemented in a minimal scheme, 
giving rise to a term sensitive to the finite parts of the hard scattering amplitudes. 

This simple symmetry argument explains the cancellations observed in [30, 31], where two- 
loop diagrams with gluon attachments to three different parton legs were shown to vanish. It 
follows that at two-loop order the soft anomalous-dimension matrix is still of the form (53). 



6.3 Three-loop analysis 

The single connected webs in Figure 6(c) containing insertions of self-energy or vertex functions 
have the same color structure as the corresponding two-loop webs and hence give nothing new 
compared with the discussion at two-loop order. This explains the cancellations observed in 
[32], where three-loop diagrams with fermionic self-energy insertions and gluon attachments 
to three different parton legs were shown to vanish. It thus suffices to consider the two four- 
gluon webs, both of which have the color structure f°-<^^f''<^'^ Tj^_ If two or more color 

generators act on the same parton (i.e., if two or more of the indices i,j,k,l coincide), then 
the products of generators with the same parton index can be decomposed into symmetric and 
anti-symmetric products. Using the Lie-algebra relations in (64), the anti-symmetric products 
can always be reduced to structures containing fewer generators. In this case we obtain one of 
the color structures already present in the two- loop case. It is therefore sufficient to consider 
only the symmetric product of four color generators, as defined in (36). We thus introduce 
the notation 

%jki = f^f" {n n Ti") + . (69) 

This object has the symmetry properties 

T^ijkl Tjilk Tikjl Tijki TkHj . (70) 
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If all four parton indices are different, then these are precisely the symmetry properties of the 
conformal cross ratio f3ijki in (50). Note that if four or three indices coincide, then the sym- 
metric product (69) vanishes. If two indices coincide, the non-vanishing index combinations 
are 

^ /f fade f bee (rparpb\ (rpcrpd\ 

J-iijj — -lijij — J J {J-i -'-%), y-i-j -'■j ). ' 

(71) 

/T- /T- /T- /T- rade rbce / rjia n~^b\ rpc rpd 

■'■iijk J-ijik jiki J-jkii / / y-l-i ■'-i ) _y_ -L j -tfc • 

The four gluons of the connected webs shown in the first two graphs in Figure 6(c) can be 
attached to up to four Wilson lines. There are thus several cases that need to be distinguished. 
If the gluons are attached to two different Wilson lines i and j, then the color structure shown 
in the first line of (71) can arise. Since two parton legs are involved in the loop diagrams, a 
dependence on fiij can arise, which is at most linear. We thus have the possibilities f3ij or a 
constant. If the four gluons are attached to three different Wilson lines i, j, and k, then the 
color structure shown in the second line of (71) can arise, which is symmetric in j and k. It 
follows that this structure can be combined with a symmetric combination of the variables , 
I3jk, and /3fc/. This leaves the three possibilities {I3ij -\- Pik), and a constant. Finally, if the 
four gluons arc attached to four different Wilson lines i, j, k, and I, then the color structure Tijki 
must be combined with a momentum structure with the same symmetry properties as those 
shown in (70). In this case four parton legs are involved in the loop integration, and hence the 
result can depend on the six cusp angles that can be formed out of the four parton momenta. 
However, the anti-symmetry of Tijki in {j, k) and (i, I) eliminates (3jk and Pu. Indeed, starting 
with any linear function of the cusp angles, symmetry arguments can be used to replace 

+ Pkl — Pik — Pjl _ ^ Pijkl 



X^ '^ijkl Pij 


— X '^ijklPkl 




{i,j,k,l) 


(i,j,k,l) 


{hj,k,l) 


%jklPik 


^ X '^ijklPjl 






{i,j,k,l) 


ii,j,k,l) 


'^ijkipjk 


= X '^ijklPil 


= 0. 


{i,j,k,l) 


{i,j,k,l) 





(i,3,k,l) 



kjki , (72) 



This leaves us with the structure TijkiPijki- 

Since conformal cross ratios are invariant under the transformation from soft to hard 

scales, the factorization constraint (49) does not prevent us from considering more complicated 
functions of such ratios. The most general possibility would be to allow a function depending 
on two linearly independent combinations of conformal cross ratios. Using (52), we can write 

F{Pijkl, Piklj — Piljk) , (73) 

where the second argument is invariant under all of the index permutations in (70). As long 
as F is odd in its first argument, F{x,y) = —F{—x,y), this ansatz respects the symmetry 
properties of the color structure Tijki- It is not easy to see how a non-trivial function of 
conformal cross ratios could arise from a calculation of Feynman diagrams. According to (12), 
any new structure in the anomalous dimension can first arise in the three-loop coefficients 
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T2 and Ti^. Since the relevant Feynman diagrams are free of subdivergences, they do not 
generate higher poles than One then expects at most single logarithms to appear in the 

coefficient of the 1/e pole term. Indeed, the structures in (72) do not contain a cusp logarithm 
and therefore can only contribute to the T2/e term in (12). An example of a structure that 
can give rise to such a term is 



fi- 



de 



+ 



-Ski 



3e 



'Sik 



3t 



-Sjl 



3t 



3 



(74) 



We thus consider a linear dependence on the conformal cross ratio Pijki as the most plausi- 
ble possibility. However, in our discussion below we will allow for an arbitrary dependence 
consistent with the symmetries of the problem. 

At this point we have exhausted the new structures that could in principle contribute to 
the soft anomalous-dimension matrix at three-loop order. Absorbing all terms not fitting the 
simple forms shown in (53) into a quantity AF^, we obtain 



(id) 



+ ^ '^iijk M^s) Pjk + f4{ois){Pij + Pik) + fbiois) 



(75) 



+ X %jkl feiois) Pijkl + F{Pijkl, Piklj — Ao'fc)] • 



{i,j,k,l) 

The functions fi{as) and F{x,y), which represents a possible dependence on conformal cross 
ratios that is more complicated than a linear term, start at three-loop order. We suppress 
the argument ctg in the latter function for brevity. Expression (75) can be simplified consid- 
erably by performing the sums over those parton indices in the second and third lines not 
involved in the various f3ij factors. In particular, the four-parton term linear in f3ijki can be 
simplified by using one of the first two relations in (72) and performing the sums over the two 
free indices. Recall that generators belonging to different partons commute, so a generator 
whose index is summed over can always be moved to the right of all other generators, and 
then relation (3) can be applied. After a straightforward calculation, we obtain 



Ar. 



X] [(/i(«s) - 2/4(q!s) - 4/6(q;J) Pij + f2{as) - 

■f %ijk fsiois) — 4/6(Q;s)j Pjk + X^ %jkl F{Pijkl, Piklj — Piljk) , 



(76) 



where the symbol ' — " means that the two expressions agree up to trivial color structures, 
which can be absorbed into (53) and hence can be dropped from ATg. Next, we evaluate the 
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constraint (49) for the above expression, which imphes 



2 \fMs) - 2h{as) - 4/6(Qi.)] + 2 '^nik [/3(«.) - 4/5(0; 



+ 



(77) 



f-i{as) - 4:fG{as) 



where in the last step wc have performed the sum over k. Relation (49) requires that this 
result be proportional to the unit matrix in color space times a coefficient depending only on 
the representation of parton i, which is not satisfied for the color structure of the first term 
in the last equation. To see this, suppose there exists a constant ki depending only on the 
representation of parton i, such that Ylj^i'^njj — 1- Taking traces over the color indices 
of either parton i, or of all other partons, we find that ki would have to be proportional to 
Cj, which evidently is not independent of the color representations of the remaining 
partons involved in the scattering process. We must therefore require that 



fsiois) = fi{as) - 2/4(^3) . 



(78) 



The most general expression for the extra terms can thus be written as (with obvious redefi- 
nitions of the coefficient functions) 



/l(as) Ai + /2(as)] + Y '^iijkfl{0is) I3jk+ ^ %jklF{Pijkl,Piklj 



iljk) 



(79) 

It follows that using arguments based on factorization and non-abelian exponentiation alone, 
one cannot exclude color and momentum structures in the soft anomalous-dimension matrix 
that are more complicated that those in (53). 

Inverting the relations between color structures that led to (76) and expressing the result 
in terms of structures containing maximal numbers of color generators, we find that the most 
general form of the additional contributions to the anomalous-dimension matrix F of n-jet 
SCET operators in (7) is 



Ar3(M,/x) 



%jkl In 



{id,k,l) 



-Sij){-Skl) 
-Sik){-Sjl) 



- f2{0is) ^ %ijk 



+ ^ %jkl F {(^ijkh f^iklj — Piljk) , 



(80) 



where the subscript "3" indicates that these structures could first arise at three-loop order. At 
this order the entire contribution AFs is proportional and all color dependence is explicit, 
i.e. it arises only from the tensors T^jki- To determine the coefficient function F{x, y) ~ 
—F{—x,y) and the numerical coefficients multiplying the other two terms it would suffice to 
calculate an arbitrary four-parton amplitude at three-loop order. If any one of the terms in 
(80) did not vanish, then our conjecture (7) for the structure of the anomalous-dimension 
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matrix would have to be modified starting at three-loop order. In Sections 6.5 and 6.6 we 
will show that the coefficients /i(as) and f2{ois) indeed vanish to all orders in perturbation 
theory, and that the function F{x, y) must vanish in all two-parton collinear limits, which is 
compatible with it being zero for all values of its arguments. We note, however, that even if the 
contribution proportional to F{x,y) would not turn out to be zero no explicit /i dependence 
enters in (80), so that 

^Ar3(W,rt = o, (81) 

and hence there would not be a contribution of this structure to the function T'{as) in (10). 
It follows that in (12) a modification could first enter in the T2/e term at three-loop order. 
Equivalently, the structure of the cusp logarithms in the anomalous-dimension matrix remains 
unaffected up to and including three loops, while the non-cusp terms remain unaffected at least 
to two- loop order. Based on our result (7), and irrespective of whether the additional terms 
in (80) vanish or not, it is therefore possible to resum large Sudakov logarithms at next-to- 
next-to-leading-logarithmic accuracy. This is sufficient for most practical applications, since 
it allows the resummation of all Sudakov logarithms which appear in next-to-leading order 
calculations of n-jet processes. Beyond this accuracy, two-loop calculations of amplitudes with 
n partons are required to obtain the necessary matching coefficients. 



6.4 Higher Casimir contributions to the cusp anomalous dimension 

For the special case of two-jet operators, the simple form (7) implies Casimir- scaling of the 
cusp anomalous dimension, i.e., the cusp anomalous dimensions of quarks and gluons are 
related to each other by the ratio of the eigenvalues Q of the quadratic Casimir operators: 

^cusp(^s) r^usp(^s) / X /OON 

— ^ = — ^ = 7cusp(as) , (82) 

see (42) and (55). This relation is indeed satisfied at three-loop order [13]. To this order 
Casimir scaling is a consequence of non-abelian exponentiation, as can be seen from our 
analysis above: restricted to the two-jet case, all possible color structures arising up to three- 
loop order are proportional to Cj. Beyond three loops non-abelian exponentiation no longer 
automatically implies Casimir scaling [9] , and there are arguments based on calculations using 
the AdS/CFT correspondence [14-16] suggesting a violation at higher orders [17-19]. The 
new color structures would involve higher Casimir invariants such as those appearing in the 
four-loop /3-function of non-abelian gauge theories [71, 72]. 

For the case of jV' = 4 SYM in the strong couphng hmit, X — g^Nc ^ oo, a, violation 
of Casimir scaling was found in [17] by considering a Wilson loop in a /c-dimensional anti- 
symmetric representation of SU (N^) in the limit where N^. and k go to infinity at fixed ratio 
Nc/k. Since the calculation was performed in the strong-coupling limit, it does not predict 
if and at which order in the weak-coupling expansion the effect would appear. On the other 
hand, it is not implausible that it might appear at some order in perturbation theory, since 
the perturbative resummation of ladder diagrams contributing to Wilson loops in J\f — A SYM 
indeed produces, when reexpanded for large A, the behavior characteristic for the strong- 
coupling limit [73]. Also, in [74, 75] an all-order form of the cusp anomalous dimension of 
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Figure 7: Four-loop connected webs involving higher Casimir invariants. 

planar A/" = 4 SYM was proposed, which is given by the solution of a certain integral equation. 
This conjecture has been checked by four-loop calculations in the weak-coupling limit [76] and 
to second order in the strong-coupling expansion using AdS / CFT and a two- loop superstring 
calculation [77]. 

Higher Casimir invariants can be constructed by considering symmetrized traces 

^aia2...a„ ^ tr [ (T^^ T^' ■ ■ ■ Tr" )+ ] 

of generators in a representation R. Any such trace contracted with n generators defines a 
Casimir invariant, since 

Cn{R, R') = d'i^"- T^^ T^' ■■■ T^' (84) 

commutes with all generators in the representation R of the group. If R is irreducible, then 
Schur's lemma implies that Cn{R,R') is proportional to the unit matrix. These Casimir 
invariants are, however, not all independent. To obtain an independent set of Casimir oper- 
ators it is sufficient to consider symmetric traces in the fundamental representation to define 
the (i-symbols, since li^^"^'""" = rf^^"^'""" with a representation-dependent index /„(/?). 

Furthermore, the invariants can be redefined, rf^^" '*" d'^'""'", such that they fulfill the or- 
thogonality conditions (f^--'^i = [78, 79]. For SU{N) groups, A^ — 1 independent 
invariants can be constructed in this way. More details on the evaluation of group-theory 
factors appearing in Feynman diagrams can be found in [80]. 

Let us now consider possible contributions of these new color structures to the cusp part 
of the soft anomalous-dimension matrix. The case n = 3 is irrelevant. The corresponding 
connected web, depicted in the middle graph in Figure 6(c), consists of three gluons attached 
to a gluon or fermion loop. These contributions have antisymmetric color structure fabc- Sym- 
metric traces of four color generators do arise, however, from the diagrams shown in Figure 7. 
The corresponding single connected webs can contribute to the soft anomalous-dimension ma- 
trix starting at four-loop order. Our goal is to study the most general contributions of these 
webs proportional to a cusp logarithm. A complete classification of potential new color and 
momentum structures that could arise at four-loop order is left for future work. 

Using the notation 

T^ijki = c^F Tf = df^'^ {T^ T/ Tf) ^ , (85) 

the possible contributions to the soft anomalous-dimension matrix linear in cusp angles have 
the following structures: f3ij "Dajj and f3ij Viuj (gluon attachments to two different Wilson 
lines), Pjkl^iijk and {Pij + (3ik)1^iijk (attachments to three different Wilson lines), or f^ijVijki 
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(83) 



(attachments to four different Wilson lines). Here we have exploited the fact that Vij^i is 
totally symmetric in its indices. Using color conservation to evaluate the sums over free 
parton indices, the result can be reduced to 

^pcusp _ j2 A, [vujj gMs) + v^j 52(«.)] + '^'^"'^ , (86) 

where the superscript "cusp" indicates that we only focus on new structures linear in cusp 
angles. The coefficient functions contain in general two terms of the form gi{as) = Uf g[ (a^) + 
gf{as), see Figure 7. They start at 0{a'^). 
Let us now evaluate the condition (49), which implies 



Q^pcusp 



= -C^{F, Ri) g2{as) + [s^w (gi{ois) - g^{ois)) + T^ijjj (g2{ois) - 2g^{as)^ 



(87) 

Only the first term on the right-hand side is of the required form and can be absorbed into 
the jet-function anomalous dimension, so that the factorization constraint (49) implies 



gi[(^s 



g2{ois) 



(88) 



The higher-Casimir cusp terms must thus have the form 



Arr" = 9i(«,) 



E 



(89) 



It is remarkable that the factorization constraint determines the structure of this term uniquely 
up to an overall coefficient function. 

The corresponding contribution to the four-loop anomalous-dimension matrix of n-jet 
SCET operators is given by 



AVr^^-giia, 



^ In ^ [Vujj + 2Vuij) + 5^ In ^ 



ijkk 



"13 



(90) 



According to (55), the cusp anomalous dimension for two-jet operators now receives a con- 
tribution not proportional to the quadratic Casimir operator of the gauge group. It is given 

by 

r:usp(«.) = Ci 7cusp(a.) - 2^i(a,) C^{F, R^) . (91) 

The four-loop cusp anomalous dimension is known for jV = 4 SYM in the planar limit [44] . 
However, we show in Appendix B that the higher Casimir contributions, as well as the three- 
loop structures (80) discussed above, are subleading in the ^ oo limit. These structures 
are therefore not visible in the planar limit. 
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6.5 A symmetry argument 

Below, we will show that (with one possible exception) the additional structures considered 
in Section 6.3 and 6.4 can be excluded, because they do not have the correct properties in the 
two-parton coUinear limit. Before turning to this discussion, however, we find it instructive to 
present a physical argument suggesting that the additional structures should be absent. To 
this end, note that the simple form of color-symmetrized soft gluon attachments to a set of 
Wilson lines discussed in Section 3.3 implies restrictions on the form of the various terms in 
(75). The Feynman rules discussed there imply that, irrespective of which Wilson lines the 
gluons attach to, one gets the same loop integral apart from substitutions of Ui vectors. 

If the symmetry properties of the diagrams also hold for the anomalous dimensions, they 
imply that the extra structures vanish due to color conservation, as we will show below. What 
makes the argument somewhat subtle is that even for IR-finite quantities, individual Feynman 
diagrams can contain IR divergences. These can manifest themselves in scaleless integrals, for 
which the expansion around d — A does not commute with the use of symmetry relations. 
For example, if one considers diagrams with exchanges between two legs i and j, one finds 
that they contain cusp logarithms, while diagrams where the gluons attach to a single leg are 
scaleless for hght-hke Wilson fines and vanish. In d dimensions, the one-leg contribution can 
be obtained by the substitutions rij — > rii and Tj — > Tj. However, after expanding around 
d = 4 the limit nj — > rii is singular for the cusp term. 

While the naive symmetry argument does not work for the cusp logarithms, we expect it to 
remain valid for those terms in the amplitude that do not depend on the momentum variables 
and fight-cone vectors of the partons in the set of Wilson lines considered in Figure 2. This 
assertion should be checked with an explicit calculation. Assuming it is true, we proceed to 
derive its implications. To this end, consider first the contributions proportional to and /g 
in the formula (75) for the most general set of extra terms arising at three-loop order. Using 
relation (72) and renaming some summation indices, these two terms can be rewritten as 

AFs 3 fsiag) ^ PijTijkk +4:f6{as) ^ I3ij%jki = h{ois) ^ XI " ^^^^ 

The last relation follows from the structure of color-symmetrized soft gluon attachments shown 
in Figure 2, when we take into account that the light-cone vectors and ni of the Wilson 
lines not involved in the cusp do not enter the value of the loop integral. It is illustrated in 
the left diagram depicted in Figure 8. Relation (92) implies that 4/6 (ct^) — fsiois)- When 
combined with relation (78), this result leads to 

fiM = Mas) - 2/4(0;,) - 4/6(a,) = . (93) 

Consider next the contributions proportional to /2 and /s in (75). They both contain two 
gluons attached to the Wilson line for parton i, plus a sum over the possible attachments of 
the remaining two gluons. In analogy with (92), we conclude that 

Ar, 3 f2{as) X Tiijj + Ma^) ^ Tujk = f2{as) ^ ^ Tujk ■ (94) 
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Figure 8: Graphical illustration of the sum over color-symmctrizcd soft gluon attachments 
giving rise to the relations /s = 4/6, gi — (left) and /2 = /s (right). 

This relation, which is illustrated in the right diagram in Figure 8, implies that 

Mas) = f2{as) - Mas) = . (95) 

We conclude that the contributions of fi{as) and M^s) to AFs in (80) vanish, which leaves 
only the possibility of a contribution involving the function F{x, y) of conformal cross ratios. 
The physical argument underlying this cancellation is the vanishing overall color charge of 
the n-parton scattering amplitude, combined with the simple form of color-symmetrized soft 
gluon attachments to collinear particles. 

For the new structure at four-loop order involving higher Casimir invariants a similar 
argument can be made. In analogy with (92), the three structures shown in (86) should all 
derive from 



2n 



2n 



mj 



(96) 



It follows that gi{as) — g2{as) — —'^gsi^s), which is incompatible with (88) unless we require 
that 

giias) = 0, (97) 

so that Ar4"'^'' in (90) vanishes. Hence, we find that the cusp anomalous dimension obeys 
Casimir scaling also at four-loop order. This observation is not in contradiction to the fact 
that, on physical grounds, one expects the finite terms in the vacuum expectation values of 
Wilson loops to receive contributions from higher Casimir invariants [9]. 



6.6 Two-parton collinear limits 

We will now rederive the conditions (93), (95), and (97) from an independent consideration. 
To this end, we analyze the behavior of the extra terms in AFs given in (80) in the two-particle 
collinear limit and check whether they are compatible with collinear factorization. For the 
contributions of the first two new structures to the anomalous dimension of the splitting 
amplitudes, we obtain 

Arsp({pi,M,/.)|^^(^^) = 2 E r,,Jin^-^^^^)('^^^^-^ 



{-Sl2){-Sij) 



+ ]nz{l - z) 



Arsp({pi,P2},)u)| 



/2(as) 



2% 



1122 



i^l,2 



(98) 
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Both contributions are incompatible with the factorization of coUinear singularities, because 
the splitting amplitudes and their anomalous dimension must not depend on the colors and 
momenta of the remaining partons not involved in the splitting process. We must therefore 
require that fi{as) = f 2(0:3) = 0, in accordance with (93) and (95). 

An analogous calculation shows that the new structure proportional to the function gi{as) 
in the expression for Ar4"^^ in (90), which would lead to a violation of Casimir scahng of the 
cusp anomalous dimension, is incompatible with the two-particle coUinear limits. Considering 
the terms proportional to ln[/x^/(— S12)] for example, we find that 



ArSp{{pi,P2},l^)\,. = -2 Vu22 + 1^1112 + Vi222 + ^^^n 



i^l,2 



In + . . . . (99) 
-S12 



The sum over T>i2n color structures cannot be expressed in terms of the color generators of 
partons 1 and 2 alone. Hence, we must require that gi{as) = 0, in agreement with (97). 

Let us then finally consider the third structure in (80). The conformal cross ratios either 
vanish or diverge when two parton momenta become coUinear. In order to study the coUinear 
limits properly, we adopt the following parameterization of the momenta of partons 1 and 2: 

2 2 

p'l ^ zEn" +p''^- , = (1 - ^)En^ -p^ - ^'^ .^ n", (100) 

4:zE 4(1 — z)E 

where — — and n • n — 2, and the ratio p±/E is a small expansion parameter. The 
coUinear limits corresponds to taking ^ at fixed E. This parameterization is such that 
Pi = P2 = remain on-shell, while — S12 = pi/[-2(l — z)]. The contribution to the anomalous 
dimension of the splitting amplitudes resulting from the last term in (80) can then be written 

as 

Arsp({pi,p2}, l^)\p= J2 FiuJij.uJij) + 4:Tuj2 F{€ij, -2(jjij)j , (101) 

where, at leading power in p±/E, 



'"^ ^'■'^ z(l — z)E \ n ■ Pi n ■ Pi j ' 

, V ^ (102) 

+ ln tfA 

Az'^{l — zYE'^ {—n-pi){—n-pj) 



Each of the two terms on the right-hand side of (101) is incompatible with the coUinear 
factorization constraint (61), unless the coefficient functions vanish in the coUinear limit pj_ 
O.'' Since the splitting amplitudes scale like l/y^si2 ^ l/|p±|) the functions must fall off at 
least as fast as e^j ~ e'^'^/^. It is an open question whether a function of transcendentality up 
to 5 (since it appears in the 1/e pole term at three- loop order) with these properties exists 
and appears in the soft anomalous-dimension matrix at three-loop order. The validity of our 
conjecture rests on the assuni])tion that F{x,y) — 0. 



^We are grateful to Lance Dixon for pointing out that a contribution to AFs in (80) involving a function 
of conformal cross ratios that vanishes in all coUinear limits in not excluded by our analysis. 
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We note in this context that a conjecture about the exponentiation of the finite terms of 
scattering ampUtudes in A/" = 4 SYM theory [76] was recently shown to be invahdated, for the 
case of n > 5 partons and at two- loop order, by Rcggc cut contributions [83], which vanish in 
all two-parton coUinear limits [84]. However, these contributions were found not to affect the 
divergent terms of the amplitude [85] . 

6.7 Extension to higher orders 

Leaving aside the possibility of functions of conformal cross ratios that vanish in all coUinear 
limits, the arguments presented in the previous sections establish our conjecture (7) at three- 
loop order and moreover exclude a certain class of modifications at four-loop order. It would 
certainly be worthwhile to test the rigor of these arguments with explicit multi-loop calcula- 
tions, as we have emphasized toward the end of Section 6.5. Nevertheless, in our opinion these 
arguments provide compelling evidence that our result is correct to all (finite) orders in per- 
turbation theory. Essentially, the constraint (49) derived from the factorization properties of 
SCET, when combined with the splitting relation (61), requires that the anomalous- dimension 
matrix must be linear in both the cusp angles and the color generators of the external par- 
tons, and that the coefficient of the cusp term is the cusp anomalous dimension. This implies 
that momentum-independent terms are color- diagonal to all orders. Momentum-dependent 
structures must have the color-dipole structure exhibited in (7). 

It thus appears that our relation (7) may indeed be an exact result of perturbative quantum 
field theory, valid in arbitrary massless gauge theories. There are few such results known in 
the literature, and it is not unreasonable to expect that the discovery of this relation will have 
profound implications for our understanding of scattering amplitudes. 

7 Summary and outlook 

We have shown that the IR poles of on-shell scattering amphtudes in massless QCD can be 
mapped onto the UV poles of the renormalization factor Z of n-jet operators in SCET. The 
RG evolution of these operators is governed by a universal anomalous-dimension matrix, whose 
form is severely constrained by soft- coUinear factorization, non-abelian exponentiation, and 
the behavior of amplitudes in coUinear limits. We have argued that only the simple form 

r(W,/x) = E ^7cusp(«.) + 5^ 

is consistent with all these constraints and have explicitly checked that they exclude any 
additional contributions up to three-loop accuracy. We also find that contributions from 
terms involving higher Casimir operators are excluded at four loops. However, our arguments 
do not exclude the presence of the term 

Ar({p}, ix) = r'^f'^' Fif^m, P^klJ - Pa,k) 
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at three-loop order (and analogous terms in higher orders), where the function F{x,y) would 
have to vanish whenever two parton momenta become coUinear, and the conformal cross ratios 
Pijki are defined in (50). We consider it unlikely that such functions arise in the anomalous- 
dimension matrix and thus conjecture that they are absent. Since the discussion in our paper 
relies solely on the commutation relations and the Jacobi identity, our results apply to any 
massless gauge theory based on a semi-simple group. Furthermore, by combining our results 
with methods developed in [86-90], which relate the singularities of massive and massless 
amplitudes, our formalism can be generalized to the massive case. This is worked out in detail 
in [91]. 

The above form of the anomalous dimension is consistent with all existing results for higher- 
order scattering amplitudes, but it would be desirable to further test it with explicit multi-loop 
calculations. It will be particularly interesting to compare with the three-loop result for the 
full four-parton amplitude in A/" = 4 SYM given in [12], once the necessary master integrals 
become available. In particular, this result will check whether color correlations between four 
partons appear, and whether they obey the constraints from soft-coUinear factorization and 
coUinear limits, i.e. whether they have the form discussed above. Also of great interest would 
be a calculation of the four-loop cusp anomalous dimensions of quarks and gluons in QCD 
or its supersymmetric extensions, either by direct calculation, extending the work of [13], 
or by using the approach based on the AdS/CFT correspondence [19]. This would test our 
prediction of Casimir scahng. The recent accomphshment of the exact evaluation of three-loop 
form factor integrals [28, 81] gives us hope that these calculations will become feasible in the 
not too distant future. 

Understanding the IR structure of scattering amplitudes is of significant theoretical inter- 
est, and having explicit results for the divergent part of the amplitudes provides an important 
check on multi-loop calculations. Also, since the singularities must cancel against those of 
diagrams with real gluon emission, our results might lead to an improved treatment of the soft 
and coUinear singularities in real emission processes. However, the most important application 
of our work are resummations of Sudakov logarithms in multi-jet processes. There is a rich 
literature on Sudakov resummation for QCD processes, starting with the pioneering papers 
[92-95] (for a review, see [96] and references therein) . In the effective theory, the resummation 
of these logarithmically-enhanced contributions is achieved by solving the RG equations for 
the Wilson coefficients. For two-jet observables, effective-theory methods have been used to 
perform resummations to N'^LL accuracy. Examples include threshold resummation for deep- 
inelastic scattering [97], Drell-Yan production [60, 98], Higgs-boson production [61, 99], and 
the extraction of from e'^e~ — > 2 jets [100]. With the anomalous dimension for the n-jet 
case at hand, it now becomes possible to reach the same accuracy also for more complicated 
observables. The evolution equation (4) is simple enough to admit exact solutions for a given 
n-parton scattering process. This can be used to perform the resummation of large Sudakov 
logarithms in closed form. 

A lot of work has been done to match parton showers with fixed-order calculations. The 
effective-theory approach allows one to not only combine fixed-order with leading-log resum- 
mations, but also to systematically resum sublcading Sudakov logarithms. To obtain predic- 
tions for n-jet observables, one needs the fixed-order results for n-parton amplitudes, which 
correspond to the Wilson coefficients of the operators in SCET. By combining our results 
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with a tree-level matrix element generator, one can obtain NLL resummations. New efficient 
methods for one-loop calculations of processes with many legs are currently used to develop 
generators for one-loop matrix elements [101]. Solving the associated RG equations then leads 
to NNLL predictions. Finally, for 2 — > 2 processes, the virtual corrections are known to 
two-loop accuracy, so that in this case N^LL accuracy can be achieved. The predictions for 
production rates of n-jet processes are obtained by combining the resummed hard-scattering 
Wilson coefficients with jet and soft functions. In contrast to fixed-order calculations or parton 
showers, these predictions are inclusive in the sense that they predict jet observables, and not 
the contributions of individual partons: the jet functions already include the integration over 
the phase space of the partons within a jet. 

For the analysis of this paper the concepts of effective field theory, such as mode separation 
and RG methods, were helpful to address an old problem of perturbative quantum field theory, 
which had resisted a solution using traditional methods. Effective field theory methods provide 
a natural language to discuss multi-scale problems, and we believe that these methods will 
play an important role in improving the accuracy of predictions for collider processes. 
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A Three-loop anomalous dimensions and Z-factor 



The three-loop expression for the renormahzation factor Z removing the IR poles of n-parton 
scattering amplitudes, as shown in (1), can be obtained by exponentiating our result (9). An 
alternative way to determine the .Z-matrix is to use the relations [102] 



2q;, — = r Z(") + (5{as) — + 



(Al) 



da.. 



91n/x ' 



where the superscript "(n)" denotes the coefficient of the l/e" pole term. These relations 
allow one to construct the higher 1/e" pole terms in a recursive way. Either way, we find 
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The expansion coefficients of the anomalous-dimensions and /3-function have been defined in 
(13). Through relations (7) and (10), the coefficients r„ and T'^ can in turn be expressed in 
terms of the expansion coefficients of the anomalous dimensions 7cusp, 7*^, and 7^, defined in 
analogy with the first relation in (13). 

We now fist these coefficients up to three-loop order in the MS renormahzation scheme. 
The expansion of the cusp anomalous dimension 7cusp to two-loop order was obtained some 
time ago [23-25, 57, 58]. The three- loop coefficient was calculated in [13]. The results are 
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The anomalous dimension 7"^ = 7'' can be determined from the three-loop expression for the 
divergent part of the on-shell quark form factor in QCD [26]. The result was extracted in [97]. 
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In the notation of this paper 27^ = 7^. We obtain 



.cl[-\^ 2.^ - 24C3) + - + 26C3) + CfTfu, (f + 1!) 

7l = Cl(~-i^^~- 68C3 + ^ a + 240(5 



^ V 4 9 135 3 ^ 3 ^ ^ 

^ V 2916 486 90 9 9 
/^2953 267r2 287r^ 512 ^ 

+ c|r,„,(^— - — - — + -C3 

^ ^ ^ / 17318 259477^ 227r^ 1928 , 

+ c.c.r.„, (^-_ + _ + _-_ (3 

^ ^2 2 /sees 40jr^ 32 , \ 

+ c.r|4(^— - — --C3j. (A4) 

Similarly, the expression for the gluon anomalous dimension can be extracted from the diver- 
gent part of the gluon form factor obtained in [26]. In terms of the anomalous dimensions 
given in [60], we have 27^'(q;s) — 7*(q;s) + ^^{ois) + l3{ois)/ois. We find 

11 4 

7o = -/^o = -y + ^ TpUf , 



/I 



„ ^2/ 692 llvr^ ^ \ /256 27i^\ ^ 

' - Ci ( -— + ^ + 2C3j +CATFnf l^—- — ] +iCFTpnf, 



g _ ^3 f 97186 6 1097r2 3197r^ 122 2071^ 
- r^29" + ^86 27r + ^^^"~9~^^~^^^^ 

^2^ /30715 11987r2 827r^ 712 
+ CAT^nf 24r + l35^^^' 

+ C.C.T.n,(fl-^-^-^C3)-2C|T.., 
^2 2 53 8 4077^ 224 \ 44 22 

Our results for 7* and 7^ are valid in conventional dimensional regularization, where polariza- 
tion vectors and spinors of all particles are treated as d-dimensional objects (so that gluons 
have (2 — 2e) helicity states). In the 't Hooft-Veltman scheme [103] or the four-dimensional 
helicity scheme [104], their values would be different. 
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B Leading-color limit 



Let us briefly discuss the iY^ — > oo limit, where a number of three- and four- loop results are 
available for jV = 4 SYM. It is interesting to ask whether these provide a test of our conjecture 
for the anomalous dimension. Unfortunately, it turns out that both the terms in (80) as well 
as the higher Casimir terms in (90) are subleading in the large-A^c limit and are thus not 
constrained by known results for planar amplitudes. 

A basis of leading color structures of n-particle gluonic amplitudes is given by the traces 
tr(i"i . . . t"") of color matrices in the fundamental representation of the gauge group [104-106] 
(see [67] for a pedagogical review). The cychcity of the trace implies that there are (n — 1)! 
different color structures for n gluons. Structures involving several traces are subleading for 
iVc — oo. These color structures can be viewed as the different possibilities of attaching n 
gluons to a quark loop. The leading contributions to squared amplitudes in the Nc ^ oo limit 
arise when a given color structure is contracted with its reverse. This gives 

TV" 

tr (r^ . . . r") tr (i"" . . . ) = ^ , (Bl) 

where terms of subleading order have been dropped. All other contractions are subleading, so 
in this sense the color basis is orthogonal. 

For Nc ^ oo, the color structure Tj • 2} acts on these traces in a particularly simple way. 
For n > 2, we find 

= -— tr(ri ... t^H^'+K.. T") , 

where again subleading terms have been dropped. Note that the result is leading order in 
color, since in the anomalous dimension this color structure gets multiplied by gi^, and g^Nc 
is held fixed in the Nc ^ oo limit. For the case of ■ Tj with i and j not adjacent we find a 
result that is of subleading order. In operator notation, we thus obtain 

Tj • Tj — ^ 6j^i±i . (B3) 
For the additional color structures arising in (80), we find 

fade f bee ^j,a j.b^^j,e j,d ^ ^^^^^ ^ 

(B4) 

jadejbce j,a j,h j,c j,d = 0{Nc), 

for all indices i,j, k, I different. Similarly, for the new structures in (90) we obtain 

Viijj = 0{N,) , Vuij = 0{N',) , Vijkk = 0{N,) . (B5) 

At n-loop order the leading color structures are of the form A^"" tr(t"i ... f""). Since the two 

structures (B4) appear first at three-loop order, while those in (B5) arise first at the level of 
four loops, their contribution is suppressed compared to the leading term. 
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